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When gauged cosmic strings form in a symmetry-breaking phase transition, the gauge field con- 
figuration at the time becomes imprinted in the spatial string distribution by the flux trapping 
mechanism. Causality and flux conservation suggest that that quantum and thermal gauge field 
fluctuations give rise to long-range superhorizon correlations in the string network. Classical field 
theory simulations in the Abelian Higgs model confirm this finding. In contrast, the Kibble-Zurek 
mechanism which most cosmic string studies are based on, only gives rise to short-distance, subhori- 
zon correlations. These results may have implications for cosmology, and it may also be possible to 
test them in superconductor experiments. 



I. INTRODUCTION 

Cosmic strings 0, [H, Q are hypothetical string-like ob- 
jects, which are stable and extremely massive and whose 
existence is predicted by many field theory [i[ and string 
theory models [5J of particle physics. If they exist in na- 
ture, they would have generally been formed in the early 
universe in symmetry breaking phase transitions or 
in brane collisions [6( in brane world models. 

Cosmic strings are cosmologically interesting, because 
their energy density tracks the overall energy density of 
the universe Most other stable objects produced 

in the early universe would eventually dominate the en- 
ergy density, in clear contradiction with observations, but 
strings always contain a small but non-negligible fraction 
of the total energy. 

Because of this scaling behaviour, cosmic strings would 
introduce nearly scale invariant density perturbations, 
and this made them at first a strong candidate for ex- 
plaining the origin of primordial density perturbations. 
Even the amplitude of the perturbations was consis- 
tent with cosmic strings in typical Grand Unified The- 
ories. However, acoustic peaks observed in the cos- 
mic microwave background [|[ and large scale struc- 
ture @ showed that perturbations existed on super- 
horizon scales. This could not be explained by the cosmic 
string model and provided the decisive piece of evidence 
in favour of inflation. 

Because cosmic strings are a fairly generic feature of 
unified theories 0], they are still being searched for in 
many different ways. If they exist, they would give a 
potentially observable contribution to the temperature 
anisotropies of the cosmic microwave background radia- 
tion [l~ol | . They would also radiate gravitational waves 
and lens light from distant objects. These effects have 
been used to constrain cosmic string models [III El, EH • 

The standard assumption in cosmic string studies is 
that the strings are formed by the Kibble-Zurek mecha- 
nism P, EH El[ • In that case, the network would be un- 
corrected at long distances, and the scaling behaviour 
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follows naturally. In numerical simulations of cosmic 
string evolution, this assumption is introduced in the 
form of random (Vachaspati-Vilenkin) initial conditions 
for the scalar field [l6j |. 

The Kibble-Zurek picture is believed to be an accurate 
description of the formation of global cosmic strings [l5| . 
However, global strings have logarithmically confining 
long-range interactions, which makes their evolution and 
cosmological effects more complicated. The usual picture 
of line-like strings with short-range interactions, which 
most cosmological studies are based on, applies only to 
gauged strings. Furthermore, cosmic strings predicted by 
realistic field theory and string theory models tend to be 
gauged rather than global [4J. 

The formation of gauged strings takes place differently 
because two fields, the scalar and the gauge field, are 
involved. Strings are formed not only by the Kibble- 
Zurek mechanism but also when gauge field fluctuations 
get trapped during the transition [17] . The predictions 
of this flux trapping mechanism have been confirmed in 
numerical simulations [TBI, El, E3| • 

This paper investigates the flux trapping mechanism 
and, in particular, the effects it has on the spatial distri- 
bution of strings. It is shown that any pre-existing fea- 
tures in the gauge field configuration become imprinted 
on the cosmic string network. If these features extend 
to super-horizon scales, which is quite possible in infla- 
tionary cosmology, they can influence the later evolution 
of the network and its observable consequences. Further- 
more, these arguments suggest that even if the gauge field 
is initially in its vacuum state, the quantum vacuum fluc- 
tuations themselves give rise to super-horizon long-range 
correlations in the string network. 



II. SYMMETRY BREAKING AND COSMIC 
STRINGS 

The simplest model of cosmic strings is the Abelian 
Higgs model in which strings correspond to Nielsen- 
Olesen vortex lines 20]. This paper will focus on this 
model, but the main conclusions are based on general 
physical principles and should be applicable much more 
broadly. They should even apply to at least certain cases 
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of cosmic superstring formation in brane inflation mod- 
els dm. 

The Abelian Higgs model consists of a complex scalar 
</> and a gauge field A^, with the Lagrangian 



:auge 



-4 V" + [D^fD^ - V(4>), (1) 



where = 9„ + ie^4 M is the covariant derivative, i 7 ^,, = 
c? M ^4^ — dyA^ is the field strength tensor, and e is the 
gauge coupling constant. We assume that the potential 
V((f>) is of the renormalisable form 



(2) 



This Lagrangian is invariant under local gauge transfor- 
mations, 

fa) _> e^W^), ^(x) - ^(s) - ^A(z), (3) 

where the rotation angle A can depend on the spacetime 
position. Physically, the field </> describes charged scalar 
particles with long-range interactions mediated by the 
gauge field. In the cosmological setting, the U(l) gauge 
group is not that of electrodynamics, although it has the 
same qualitative properties. In this paper, the compo- 
nents of the field strength tensor will nevertheless be 
referred to as "electric" and "magnetic" fields. 

In the limit e — * 0, the gauge field decouples, and the 
Lagrangian reduces to 



£ g iobai = d^*d»4>-v{4>). 



(4) 



This Lagrangian is only invariant under global phase ro- 
tations (f> — ► e lA (j), where A is constant over the whole 
spacetime. Physically, the quanta of </> are therefore 
"electrically" neutral. 

In both cases, when m 2 is positive, the vacuum state, 
given by the minimum of the potential, corresponds to 
= 0. For negative m 2 , there is a set of degenerate vacua 
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m 
2A 



(5) 



parameterized by the complex phase — tt < 8 < it. Be- 
cause these vacua are related by a symmetry, they are 
all identical. Once the system has chosen a particular 
vacuum, its state is no longer invariant under symmetry 
transformations, and the symmetry is said to be spon- 
taneously broken. In the gauge theory, the covariant 
derivative term gives the gauge field a mass m 7 = ev. 

In both models there are line-like stable topological 
defect solutions in the broken phase. The solution has 
the form (in cylindrical polar coordinates r, if, z) 



(r,ip,z)=vf(r)e iN ^, 



(6) 



where N w is an integer winding number, and f(r) is a 
function that has to be determined numerically and sat- 
isfies the boundary conditions 



f(r) 



0, as r — > 0, 

1, as r — » 00. 



(7) 



In the gauge theory, the gauge field in the appropriate 
gauge is 



A 



N„ 



■a{r)ip, 



(8) 



where a(r) is a function that has to be determined nu- 
merically and satisfies the same boundary conditions as 
f(r). Because the magnetic flux is given by the contour 
integral 



r 2tt 
<T> = •'/> df-A = — N u 



(9) 



the strings carry a quantized magnetic flux. In fact they 
are fully analogous to Abrikosov flux tubes in supercon- 
ductors [2ll| . Outside the vortices, the magnetic field 
decays exponentially, as a manifestation of the Meissner 
effect. 

In the gauge theory |T]), the energy of the string is 
localised exponentially on the string, and these gauged 
strings behave like idealised Nambu-Goto strings [2. |22|. 
In contrast, energy of a string is not localised in the global 
theory ([4]). The strings have a confining logarithmic in- 
teraction, and therefore these global strings evolve differ- 
ently, and much less is known about their cosmological 
effects. 

If a network of cosmic strings was formed in the early 
universe, its later evolution will have produced poten- 
tially observable effects such as gravitational wave emis- 
sion and temperature anisotropics in the cosmic mi- 
crowave background (CMB) radiation. All the observ- 
able effects are gravitational and depend therefore on the 
dimensionless combination Gfi, where G is Newton's con- 
stant and fi is the string tension, defined as the energy 
per unit length. 

The evolution of the string network is a complicated 
non-linear problem, and even with large-scale numerical 
simulations one can follow it only for a relatively short 
time [TTj . To make predictions, one usually assumes that 
the evolution is such that the network scales with the 
horizon size [23|, [24[ ■ As will be discussed in more detail 
in Section IVlIl this scaling behaviour follows naturally if 
there is sufficient energy loss and if the string network is 
uncorrelated at long distances. 

The scaling behaviour means that at any time, the net- 
work looks statistically the same in units of 1/H. With 
this assumption, the total length of the string network 
inside one Hubble volume scales as H^ 1 , indicating that 
the total energy is ~ fi/ 1 H. Dividing by the volume 
~ H~ 3 , we find the energy density, 



-Gftpc, 



(10) 



where p c is the critical energy density. Therefore as long 
as the scaling assumption is valid, the strings contribute a 
constant fraction of the total energy of the universe, and 
will never dominate or be negligible. This is why they 
are cosmologically more interesting than other topologi- 
cal defects, i.e., domain walls or monopoles, which tend 
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to dominate the energy density and overdose the uni- 
verse. 

As a consequence of the scaling, the perturbations gen- 
erated by the string network are scale invariant as re- 
quired by the observations. However, in the absence of 
super-horizon correlations, they would not have the ob- 
served acoustic peaks (25|. The contribution to the den- 
sity perturbations from cosmic string must therefore be 
subdominant. This sets an upper limit for the string ten- 
sion Gp < 10~ 5 . 



III. LOCAL STRING DENSITY 

To discuss correlations in the string network, Liu and 
Mazenko [26[ defined the local string density p as 



p(x)=^2 I ' ds-^5(x- x a (s)) . 



(11) 



where a labels the strings and s is a coordinate along the 
string. In the Abelian Higgs model, this can be expressed 
in terms of the fields as 



P = V x q, 



where q is defined as 



— ( Im 6*D(i) - el) = — V0. 



(12) 



(13) 



Here 8 = arg0 is the complex phase of the scalar field 
(f>, and (j> — (j>/\4>\ — exp(i6). The local string density is 
gauge invariant and sourceless 



V-p = 0. 



(14) 



It vanishes outside the strings and has a delta-function 
peak tangential to the string at the centre of the string. 

The two-point string correlation function Gij (x) is de- 
fined as 



G ij (x) = (p i (0)p j (x)). 



It is also often useful to consider the Fourier transform 
of this, 

Gij{k)= f d 3 xe lls G tJ (x), (16) 



which gives the momentum space correlation function 

(pi(k) Pj {k')) = G ij (fe)(2 7 r) 3 <5 (* + *') . (17) 

Because of rotation invariance and the sourcelessness 
condition (fT4|) , this can be written as 

G ij (k) = U j - 1 ^\G(k), (18) 



where G(k) is a scalar and can be obtained from the trace 
of G M , 



G(k) = -G u (k). 



(19) 



Note that the correlator G is not the same as the corre- 
lator of the energy-momentum tensor [IJ, [2?], [28| , which 
is more commonly discussed and which seeds density per- 
turbations. Long-range correlations in one quantity do 
not imply long-range correlations in the other. However, 
as will be discussed in Section IVIII long-range correla- 
tions in p will affect the time evolution of the network 
and thereby also indirectly the density perturbations and 
other observable signatures. 



IV. STRING FORMATION 

Cosmic strings were formed in the early universe if a 
U(l) symmetry became spontaneously broken in a phase 
transition. In the original scenario [l|, the transition took 
place at a high temperature. To first approximation, 
the effects of the temperature are given by an additive 
temperature-dependent contribution to the mass param- 
eter m 2 , so that the effective mass parameter m 2 (T) is 



m 2 (T) 



m 



A T 



where m 2 < 0. At high temperatures, 



T > Tc 



3e 2 + 4A' 



(20) 



(21) 



the effective mass parameter is positive and the symme- 
try is restored. When the universe cooled down and the 
temperature decreased below the critical value, m 2 (T) 
became negative and the symmetry was broken. 

In inflationary models, string formation can take place 
at the end of inflation, and in this case the temperature is 
generally extremely low. In hybrid inflation [29(, the field 
(15) is coupled to the inflaton field a, so that the effective 



2/ \ 2,22 

m [a) = m + g a , 



(22) 



where g is the coupling constant for the <f> — a interac- 
tion. When the inflaton rolls below the critical value 
u c = {— m 2 / g 2 ) 1 / 2 , the symmetry gets broken. This tran- 
sition not only forms strings but also triggers the end of 
inflation. 

Brane collisions in brane world and string theory mod- 
els are effectively very similar to the phase transition in 
hybrid inflation [§[. In general, the role of the inflaton is 
played by the distance between the branes, and when the 
branes collide, they annihilate through tachyon conden- 
sation, which is effectively a symmetry-breaking phase 
transition. 
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Whatever the microphysics that causes the symmetry 
breaking, it generally forms cosmic strings. There are 
two distinct physical mechanisms that are responsible for 
this: the Kibble-Zurek mechanism [J [H| and flux trap- 
ping [II, [13. 



A. Kibble-Zurek Mechanism 

The Kibble-Zurek mechanism involves only the scalar 
field 0, and therefore it operates in both global and gauge 
theories. In the gauge theory, it is only part of the dy- 
namics and becomes dominant if the gauge field can be 
ignored. 

When the symmetry is broken, the scalar field has to 
choose one of the possible vacua (O, corresponding to one 
particular complex phase 9. All values of are equally 
likely because all the vacua are identical. Kibble's origi- 
nal observation jj was that in cosmology, any two points 
that are separated by more than the horizon size are 
causally disconnected, and they will have to make this 
choice independently of each other. After the transition, 
the universe will therefore consist of finite regions, which 
cannot be larger than the horizon, in each of which the 
choice of the vacuum is random, and between which it is 
uncorrelated. 

If one follows the complex phase around a closed curve 
that passes through at least three of these regions, there 
is a certain non-zero probability that it changes by a 
multiple of 2tt. In this case, continuity implies that the 
field has to vanish somewhere inside the curve, and this 
corresponds to the core of the string. According to this 
argument there should be roughly one string per corre- 
lated region, whose size is limited by the horizon size. 

Zurek [l4[ showed later that the size of the correlated 
regions is generally determined not by the cosmological 
horizon size, but by the critical slowing down during the 
transition. The correlation length £ grows as the critical 
temperature is approached, and in an adiabatic transi- 
tion it would diverge at the transition point. However, 
the phase transition can never be adiabatic. The dynam- 
ics of the system slows down when the transition is ap- 
proached, and if the transition takes place in a finite time, 
the system cannot stay in equilibrium. Instead of actu- 
ally diverging, the correlation length remains finite, freez- 
ing out at some finite value £, which determines the size 
of the correlated regions. The ultimately upper bound 
for this is given by causality, since the correlation length 
cannot grow faster than the speed of light. From this 
argument Zurek was able to derive a scaling law for the 
number of defects as function of the cooling rate, which 
has been confirmed in experiments [l5l |. 

Since the complex phase 9 is only correlated at dis- 
tances less than £, Eqs. (fT2")) and (jT3J) imply that the 
same must be true for the string density p. Therefore 
the correlator must have a finite range. For instance, Liu 



and Mazenko find a Gaussian correlator [26j 

G tJ (x)^e- s2 /^ 2 . (23) 

Whatever the precise form of the correlator, there are no 
long-range correlations in a string network formed by the 
Kibble-Zurek mechanism. 

Furthermore, the integral of the correlator over the 
whole space vanishes. It can be expressed as a surface 
integral over its boundary, 



d^xGijix) = (piiO) / d d xpj(x) 



e jk i I d 2 S k (p t (0) qi (x)). (24) 



In this integral all points x are on the boundary, and 
because of the finite correlation length, the correlator 
(pi(Q)qi(x)) and therefore the right-hand side of Eq. (|24|) 
vanishes, i.e., 



d 3 xG i3 (x) = 0. 



(25) 



This indicates a negative correlation between strings: 
The winding number of a string is screened by the wind- 
ings of nearby strings. 

Because of the finite range, the momentum space corre- 
lator G(k) is analytic in k 2 and can be Taylor expanded, 



n poo 

G{k) = — I drrsm(kr)G u (r) 



10 

A f°° , 2 / k 2 r 2 

47r / drr\l 

Jo V 6 

~2 , nnA 



G{r) + 0(fc 4 ) 



£fc^ + 0(fc 4 ). 



(26) 



The constant term cancels because of Eq. (f25j) , and the k 2 
term has to have this form because £ is the only available 
dimensionful scale. 



B. Flux Trapping 



ping 



In the gauge theory, strings can also form by flux trap- 
Magnetic fields cannot exist in the bro- 



lf 



17] 



ken phase vacuum, and if a non-zero field is present at 
the time of the transition, it must therefore become con- 
fined in strings. The simplest example is a transition in 
a uniform magnetic field, which is a common setup for 
superconductor experiments. The Meissner effect tries to 
expel the field but inside the superconductor it is unable 
to do that, and instead it forms a lattice of Abrikosov 
vortices [21] , where the average density of strings is di- 
rectly related to the field strength, 



(27) 



In the early universe there was no external field, but 
there were thermal and quantum fluctuations. Because 
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the photon is massless in the symmetric phase, these fluc- 
tuations can have arbitrarily long wavelengths, or equiva- 
lently, arbitrarily low wave numbers k. When the system 
enters the broken phase, magnetic field modes with wave 
number less than the photon mass, |fc| < m 7 , are pro- 
hibited. The system tries to suppress these modes to 
minimize its energy, but the magnetic flux is conserved 
as a consequence of the Maxwell equation 

— + V x E = 0. (28) 
ot 

This means that the net flux through a surface can only 
change by a boundary contribution, when magnetic field 
lines move across the boundary, 

— = -jdr f -E. (29) 

Because of this, modes with longer wavelengths take more 
time to decay. Essentially a mode with wavenumber k 
cannot decay in shorter time than what it takes for light 
to travel one wavelength A = l/\k\. 

A long-wavelength magnetic field mode that has not 
had time to decay before the phase transition is com- 
pleted remains in the system even though it is now in the 
broken phase. As with a uniform field discussed above, 
this is only possible if it is confined inside strings. The 
string density (averaged over suitable volume) is then lo- 
cally related to the magnetic field strength by the same 
equation (j2"T|) . 

This means that long- wavelength fluctuations break up 
into strings in the same way as a uniform field. This 
is natural because a fluctuation with long enough wave- 
length (for instance longer than the cosmological horizon 
size) is locally indistinguishable from a uniform field. 

On the other hand, short-wavelength fluctuations be- 
have adiabatically and decay, which means that there 
must be a critical wavenumber k separating the two dif- 
ferent behaviours: Modes with |fe| > k decay but modes 
with \k\ < k become trapped into strings. The critical 
wavenumber depends on the dynamics of the system and 
is generally lower for slower transitions. In principle, its 
precise value can be calculated for any particular sce- 
nario, but that will not be necessary for the conclusions 
drawn this paper. 

If there was any feature in the initial magnetic field 
configuration that was larger than the critical length 
scale 1/fe, this feature would survive the transition and 
be imprinted in the cosmic string network. For instance, 
a region with a strong uniform magnetic field would cre- 
ate a bunch of nearly parallel cosmic strings, which would 
roughly follow the initial magnetic field lines. Even in the 
absence of such specific features, which are unlikely in 
typical cosmological models, thermal and quantum fluc- 
tuations lead to non-trivial effects, and as discussed in 
the next section. 



V. CORRELATIONS FROM FLUX TRAPPING 

Because the magnetic field is sourceless, its momentum 
space two-point correlation function has the same tensor 
structure as Eq. (fTS)) , 

(BiikjBjik')) = G? 3 (k)(2ir) 3 6 (k + fc') , (30) 

where 

Gf J (k)=(s i3 -^G B (k), (31) 

and G B (k) is a scalar function of k = \k\. The coordinate 
space correlator G B (x) is given by the Fourier transform 

GfjW = J -^^ s G B {k). (32) 

The conservation of magnetic flux (|2"9")l means that the 
correlator can only change by a spatial derivative. In 
Fourier space, the change in G B is at least quadratic in 

k, 

G B (t; k) = G B (0; k) + (fc 2 <% - fe i fc i )A(fc). (33) 

Once the system enters the broken phase, the long- 
wavelength magnetic field modes freeze out as discussed 
in Section HVBI The final correlation function G B nal will 
therefore interpolate between the initial on G 3 ^ at low 
wave number k < k and the equilibrium correlator G B q 

at high wave number k > k. 

The precise form of the correlator is not important, but 
one obtains a useful approximation by assuming that the 
long-distance dynamics are diffusive. Then each Fourier 
mode decays as 

B(t,k) = B(0,k)e- k2t/47ra 1 (34) 

where a is the electric conductivity. At some time t the 
dynamics become non-linear when strings form, and one 
can identify 

* - (35) 

so that the final correlator is 

Gf nal (fc) = Gf ni (fc)e- fc2 /2 p (36) 

When the dynamics become non-linear, these magnetic 
field fluctuations are turned into cosmic strings. If one 
ignores the contribution from the Kibble-Zurek mecha- 
nism, so that flux trapping is the only source of strings, 
Eq. (|27f implies that the string density correlator is sim- 
ply given by 

G ( fc ) = T^fnalW = ^Gf ni (fcK fc2 / 2p . (37) 
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To calculate the correlations of strings formed by flux 
trapping, one needs to know the spectrum of initial fluc- 
tuations G^;, which depends on the scenario. The two 
simplest and most relevant cases, when the field in ini- 
tially in vacuum or in thermal equilibrium, are discussed 
below. 



A. Vacuum Fluctuations 

The simplest situation is when the gauge field is in 
its vacuum state at the time of string formation. This 
is the case after a long period of inflation, which would 
have wiped out any excitations. Because the gauge field 
is conformally invariant, no long-range fluctuations are 
generated during inflation [3(J. Therefore the expansion 
of the universe simply redshifts and dilutes any existing 
gauge field configuration, taking it to its vacuum state. 

In this case, the symmetry breaking phase transition 
can be thought of as a non-adiabatic change of the vac- 
uum state of the system. Normally this leads to pair 
production of particles, which is often be described using 
Bogoliubov transformations within linear theory with a 
time-dependent Hamiltonian. Our case is physically very 
similar but the produced excitations are cosmic strings, 
and because they are non-linear objects, the linear ap- 
proximation is not valid. In contrast, the physical pic- 
ture presented in Section [Vj which is based on causality 
and conservation laws, should still be applicable. 

It is important to note that even if the universe re- 
heated between the end of inflation and the phase tran- 
sition, the cosmologically relevant modes would still re- 
main at zero temperature, because they are outside the 
horizon and can therefore not be affected by local reheat- 
ing dynamics. Therefore the calculations in this Section 
should also apply to that situations. Furthermore, these 
calculations are not restricted to cosmological scenarios 
but apply generally to any system at zero temperature. 

In the symmetric phase, the vacuum two-point corre- 
lation function of the magnetic field can be calculated 
from the photon Feynman propagator, 



IT} 



k 2 



(38) 



where we have adopted the Feynman gauge. This is unaf- 
fected by the expansion of the universe because the gauge 
field is conformally invariant. The equal-time two-point 
correlator of the gauge field is 



2tt 



S. 



(39) 



The magnetic field correlator is obtained by calculating 
the curl of the gauge field, Bi(k) — ieijkkjAk(k), giving 



G B {k) 



k 
2' 



(40) 



The coordinate space correlator is given by the Fourier 
transform of this, and has a power-law form 



G B (r) = - : 



1 



(41) 



This means that the quantum vacuum fluctuations of 
the magnetic field are correlated over infinitely long dis- 
tances. 

It is important to understand that the non-zero space- 
like two-point function (|4ip does not indicate superlu- 
minal signal propagation and violate causality. It de- 
scribes how the field fluctuations are correlated, not how 
they propagate. In cosmology, one would normally ar- 
gue that nothing can be correlated over more than the 
particle horizon. However, the way inflation solves the 
horizon problem is precisely by making the particle hori- 
zon larger than the currently observable universe. The 
power-law form (|4"Tj) should therefore be valid at all ob- 
servable scales. 

In the broken phase, the Higgs mechanism makes the 
field massive and the vacuum two-point function is 



G B {k) 



2jk 2 



(42) 



where m 7 is the photon mass. The coordinate space cor- 
relation function can be obtained by taking the Fourier 
transform of this. 



G B (r) 



4 7r 2 r 2 



m-,r + 



(43) 

where K v is the modified Bessel function of the second 
kind. The asymptotic long-distance behaviour is expo- 
nentially decaying, 



2 l/2 



5/2 



8?r 3/2 r 3/2 



(44) 



In an adiabatic transition, the correlator would change 
from Eq. (|40f to Eq. I]42p. and the long-range correla- 
tions would disappear. However, as was discussed in Sec- 
tion IIVB1 this cannot happen in a finite time. Instead, 
flux conservation suggests that the longest wavelengths 
must freeze out and retain their original amplitude (|40|) . 
This is much higher than the amplitude of vacuum fluc- 
tuations in the broken phase (|42p , and therefore the tran- 
sition has turned the quantum vacuum fluctuation into a 
classical magnetic field. 

According to Eq. ([56]) . the frozen-out magnetic field 
correlator after the transition is 



final 



(k) 



-k 2 /2k 2 



(45) 



where the broken phase vacuum contribution (|42p has 
been dropped because to 7 is generally larger than the 
scales that are relevant here. By causality, the critical 
wavenumber k has to be greater than the Hubble rate 
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|fc| > H. The coordinate space correlator is obtained by 
Fourier transforming this, 



final 



(r) = 



4tt 2 



1- 



D -r 2 k 2 /2 



— ( — — rk | erfi [ — = . 
2 Vrfc J \V2j 

(46) 

where erfi is the imaginary error integral. Unsurpris- 
ingly, since the long-distance modes are unchanged, the 
long-distance asymptotic behaviour remains the same as 
before the transition 

1 



2ir 2 r 4 



(47) 



According to Eq. (|3T() , these correlations are imprinted 
in the string network. Including the contribution from 
the Kibble-Zurek mechanism (f2"6")) . Eq. (|4"5)) predicts 
that the string density correlator should behave at low 
wavenumbers as 



£fc 2 + 0(fc 3 



(48) 



Similarly, Eq. (|47| implies that the asymptotic long- 
range correlator in coordinate space is 



G(r) ~ - 



8?r 4 r 4 ' 



(49) 



The strings are therefore correlated even on superhorizon 
scales, even if the initial state was empty vacuum. 



B. Thermal Fluctuations 

The correlations are even stronger if the long- 
wavelength modes of the gauge field were initially not 
in vacuum. Reheating or any other post-inflationary 
phenomena cannot create these conditions, because by 
causality they can only affect sub-horizon modes. How- 
ever, if inflation lasted only the 50 or 60 e-foldings re- 
quired to solve the horizon and flatness problems, pre- 
inflationary field fluctuations would have been imprinted 
in the cosmic strings. In some models, trans-Planckian 
physics or phenomena taking place during inflation can 
also take the gauge field away from its vacuum state. 

For concreteness, let us assume that before the start 
of inflation, sub-horizon magnetic field modes were in 
thermal equilibrium at temperature Tj n ;. During infla- 
tion, these comoving modes left the horizon and their 
temperature was redshifted down to T = e~ N T^, where 
N w 50 is the number of e-foldings. Today, the wave- 
length of these modes is comparable to the horizon size, 
and they are therefore be in the observable range. 

At the time of string formation, at the end of inflation, 
the relevant long-wavelength field modes were therefore 
in thermal equilibrium at temperature T. The equal-time 
finite-temperature gauge field correlator in the Feynman 
gauge is 



(A t (k)A,(k')) = 6^(2^ 5(k + k>) A coth A. 



(50) 



This gives the magnetic field correlator in momentum 
space, 



k 



k 



G*(fc) = -coth— , 



2 2T 

and, through a Fourier transform, in coordinate space 



(51) 



x 



T 2 ( , 2 cothTrTr 

x — - 1 - coth 2 ixTr 

4r 2 \ nTr 



c- X{Xj 

% - 1 x 



ttT 3 
"27 



• coth irTr ( 1 - coth 2 vrTr ) 



(52) 



where r = \x\. The long-distance behaviour is obtained 
by taking r — > oo, 



G\ 



|(f) -(i 



! 1 J 



Airr 3 J 



(53) 



Taking the low wavenumber limit of Eq. (|5ip and using 
Eq. (I37|) . one find the prediction for the string density 
correlator, 



G(fc)«^A + ^ 2 + 0(fc 4 ). 



(54) 



In coordinate space, Eq. (fBU)) implies that the long- 
distance correlator is 



G 



XiXj 



e 2 T 
16n 3 r 3 



(55) 



This is traceless, indicating that G(r) — Gu/2 decays 
exponentially and is therefore not a particularly useful 
quantity for discussing the correlations. Nevertheless, 
even the trace behaves very differently from the Kibble- 
Zurek case because winding number is not screened, 



d 3 xG(x) = lim G(k) 



e 2 T 



i ... 4ir 2 



(56) 



which should be compared with Eq. (j25|) . 

Comparison of Eq. (|55[) with Eq. (|4"9")l shows that the 
correlator decays more slowly than at zero temperature. 
Because the cosmologically relevant distance scales are 
extremely long, this means that even a very low non-zero 
temperature can have significant effects. 



VI. SIMULATIONS 

Predictions of the the Kibble-Zurek mechanism and 
flux trapping have been tested and confirmed in numer- 
ical simulations [H [H 0J, [M H . Unfortunately these 
tests are limited to classical field theory, because there 



are no good ways of simulating non-equilibrium dynam- 
ics in quantum field theories, except in the simplest mod- 
els [33]. For example, the two-particle irreducible ef- 
fective action formalism (34|, which is a promising tool 
for perturbative non-equilibrium processes, fails to incor- 
porate topological defects [35[. Interestingly, it is pos- 
sible to test these effects experimentally in condensed 
matter systems which as superfluids and superconduc- 
tors m m, m, m m sn, m, s m, which are, of 

course, fully quantum mechanical. So far, the emphasis 
has been on determining the number of strings, and the 
tests of spatial correlations have so far been limited to 
short-distance effects [H, [H, ||U . 

To test the predictions in Section |V] I carried out a 
set of classical field theory simulations of the gauge the- 
ory (fT]). The simulations are perfectly conventional, de- 
scribing the field evolution of the Abelian Higgs model 
in expanding space. The only difference from usual cos- 
mic string simulations are the initial conditions, which 
were calculated from the actual thermal equilibrium or 
vacuum state of the quantum field theory, rather than 
imposing the usual ad-hoc random phase initial condi- 
tions Hal. 
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FIG. 1: The string density correlator G(k) measured at the 
end of the simulation with thermal initial conditions. The 
empty and filled diamonds correspond the data for e = 0.5 
and e = 0.1, respectively, and the circles are the correspond- 
ing data after cooling by gradient flow. The two solid lines 
show the flux trapping prediction (|58[) for e = 0.5 (upper) 
and e = 0.1 (lower). The dashed and dotted lines include an 
additive contribution £fc 2 , as predicted by the Kibble-Zurek 
mechanism (|54l) . 



A. Non-Zero Temperature 

The first set of simulations focused on the thermal case 
discussed in Section fVBl They were carried out on a 256 3 
lattice with comoving lattice spacing 5x = 1. The scalar 
coupling constant was A = 1, and for the gauge coupling 
two different values e = 0.1 and e = 0.5 were used. 

Initially, the system was prepared in classical thermal 
equilibrium at temperature T = 0.5, corresponding to 
the canonical ensemble p cx exp(—/3H). This was done 
with a Monte Carlo algorithm, more details of which can 
be found in Ref. [47j . 

The fields were then evolved according to the classical 
equations of motion. This was done in conformal time 
r defined by dr — dt/a(t), where a is the scale factor. 
In terms of rescaled fields <j> = a<\> and Ei = —d T Ai, the 
equations of motion are 



d 2 T 4> 
d r Ei 
diEi 



= D 



{m 2 a 2 - d 2 T a/a)4>-2\\(t>\ 2 <i>, 
2e 

2eIm</>*9 T <; 



djFij + 2< Una' I) 



(57) 



The universe was assumed to be radiation dominated, so 
that a = 1 + i?iniT, where = 0.1 is the initial Hubble 
rate, and consequently d 2 a/a = 0. The only direct effect 
of the expansion is therefore the growth of the mass term. 
These equations were discretised with constant conformal 
time step St = 0.05. 

To allow a meaningful comparison between the two 
values of e, the mass parameter m 2 was set to m 2 = 
-0.285 for e = 0.5 and to m 2 = -0.225 for e = 0.1, so 
that in both cases the transition to the broken phase took 
place near a f=a a/2- 



The local string density was measured at scale factor 
a = 2 using the gauge-invariant lattice winding num- 
ber (48[, and the two-point correlation function G(k) was 
calculated. The results are shown as diamonds in Fig. Q] 

Because at a = 2 the temperature is still fairly high, 
T = 0.25, thermal fluctuations are present especially at 
short distance. They should not affect the long-distance 
behaviour, but to confirm this, the circles in Fig. [1] show 
the same data after cooling the system by gradient flow. 
This removes the fluctuations and makes sure that the 
strings can really be thought of as classical solutions. 

The solid lines show the predicted contribution from 
flux trapping, 



(58) 



which seems to describe the asymptotic behaviour cor- 
rectly. 

The dashed and dotted lines show the predicted long- 
distance behaviour (|5~i|) , with Kibble-Zurek correlation 
length £ determined as the best fit to the data at k < 0.1. 
Without cooling, its value was 



0.226 ±0.016, fore = 0.5, 
0.127 ±0.003, fore = 0.1, 



and after gradient flow 



0.61 ±0.01, fore = 0.5, 
0.66 ±0.01, fore = 0.1. 



(59) 



(60) 



As one would expect, the correlation length grows during 
cooling, as the field becomes more ordered, but the long- 
distance behaviour is unchanged. Note that having £ less 
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than one does not pose a problem for the lattice discreti- 
sation, because £ is not the actual correlation length but 
only a fit parameter that is proportional to it. 

The agreement of these predictions with the data at 
long distances (low k) is remarkable. This supports the 
calculations in Section [V] and suggests that the contri- 
butions from the two mechanisms are additive to a good 
approximation. One can therefore consider the two mech- 
anisms to be essentially independent of each other. 

In Fig. [1] there appears to be a discrepancy between 
the prediction ([54|) and the data at short distances, as 
the measured correlator turns down at high k. However, 
this does not indicate any conflict between the theory 
and the data because Eq. (|54|) is only the asymptotic 
long-range form, and a drop at high k is to be expected. 

B. Vacuum 

The second set of simulations tested the zero- 
temperature case discussed in Section IV Al Following 
an approach that is widely used in studies of inflation- 
ary preheating [49|, [5(J HH, quantum fluctuations were 
modelled with an Gaussian ensemble of classical fluctua- 
tions with the same two-point function. This should be 
a reasonable approximation to model quantum dynamics 
by classical fluctuations with the same initial two-point 
correlation function (49|: At linear level, which should 
be a reasonably good approximation at early times, the 
classical and quantum dynamics are identical, and at late 
times, the occupation numbers of the relevant modes are 
high and therefore classical approximation should again 
be valid. 

Because the predicted correlations are weaker in the 
vacuum case than at non-zero temperature, larger lat- 
tices are needed. To save memory and computing power, 
and elongated lattice of size 128 x 128 x 32768 was used 
and the correlations were measured in the long direction. 
Because even the shorter directions are longer than the 
light-crossing time, the shape of the lattice should not 
affect the measured correlations in any way. 

Again, the field evolution is given by Eq. (|57p . but 
because the initial state is in vacuum, the expansion alone 
would not cause a transition in the full quantum theory. 
It is therefore necessary to have some mechanism that 
changes the effective mass term, such as coupling to the 
inflaton field in hybrid inflation. For simplicity, the scale 
factor was set to a constant, a = 1, and the mass term 
had linear time dependence 

m 2 (r) — mf ni — ct. (61) 

The mass term was initially set to mf ni — 1, and it was 
quenched at rate c = 0.2 for time At = 10, so that the 
final mass term was ?ri^ nal = — 1. 

The string density correlator measured immediately af- 
ter this quench is shown in Fig. UK a), and the same mea- 
surement after gradient flow cooling in Fig.^b). Because 
of the anisotropic lattice shape, the correlator was only 



(a) 




FIG. 2: The string density correlator G(k) measured with 
vacuum initial conditions (a) at the end of the quench and 
(b) after gradient flow cooling. The empty and filled circles 
correspond the data for e = 0.4 and e = 0.1, respectively. 
The two solid lines show the asymptotic behaviour predicted 
by flux trapping (I63|l for e = 0.4 (upper) and e = 0.1 (lower). 
The two dashed lines include an additive contribution £k 2 , as 
predicted by the Kibble-Zurek mechanism (148 [) . 

measured in the long (z) direction, i.e., the measured 
quantity was 

G(k) = ^(G xx {kz)+G vy (kz)). (62) 

The solid lines show the asymptotic prediction from 
flux trapping only, 

G <Mi)1< < 63 > 

and again they agree with the measured asymptotic be- 
haviour. 

The dashed lines show the predicted long-wavelength 
correlator (|48p. which includes contributions from both 
flux trapping and the Kibble-Zurek mechanism. Imme- 
diately after the quench [Fig.^a)], the best fit values of 
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the Kibble-Zurek fit parameter £ were 

- f 0.154 ± 0.002, fore = 0.4, 
jo.l31± 0.003, fore = 0.1, 

and after gradient flow [Fig. E^b)], 



(64) 



0.812 ±0.010, 
0.605 ± 0.020, 



for e = 0.4, 
for e = 0.1. 



(65) 



As in the thermal case, the Kibble-Zurek correlation 
length grows during gradient flow as the system gets more 
ordered and the smallest string loops disappear, but the 
long-range asymptotic behaviour remains unchanged and 
is correctly predicted by the flux trapping mechanism. 



VII. COSMOLOGY 

Most calculations of cosmological effects of cosmic 
strings rely on scaling. The strings were formed in the 
very early universe and one needs to know their whole 
cosmological evolution. Because of the vast time and 
length scales involved, it is impossible to calculate this in 
detail. However, it is believed that the statistical proper- 
ties of the string network scale in such a way that, relative 
to the horizon size, they look the same at all times Q. 

By considering the string density correlator G(k), one 
can see that this is natural if the strings were formed by 
the Kibble-Zurek mechanism. At low k, the correlator 
grows as k 2 , and at some scale it turns down and falls as 
1/k. If we label the position of the maximum by fc max , 
we can approximate the correlator by 



G(k) 




for k < k max , 
for k > fc maY . 



(66) 



Expressing this in units of H, one has 



l,/ \(HO(k max /Hf/ K , 



for k < k n 

for K > kn 



c/H, 
</H, 



(67) 

where k = k/H. 

Because of causality, the correlation length £ cannot 
grow faster than the horizon size l/H. Assuming that it 
saturates this bound, the first coefficient Ht; is constant. 

As the string network evolves, its loses energy by ra- 
diating it into other degrees of freedom. Because of this, 
string loops shrink and eventually annihilate. This sup- 
pressed the string density correlator, but because this 
process can only take place on subhorizon scales, it only 
effects high wave numbers k > H. Again, if one assumes 
that this energy loss is efficient so that no power accu- 
mulates at high k, one concludes that the peak of the 
correlator moves as fc max cx H. In Hubble units, the 
correlator (167|) would therefore remain unchanged, which 
means that the network scales. 



The introduction of the thermal flux trapping contri- 
bution clearly violates this scaling. Eq. (|6"Tj) becomes 



G(K)/H: 



\ + (#£V, for k < k max /H, 

\ (Hi)(k max /H) 3 /K, for k > k max /H, 

(68) 

where T s and a s are the temperature and the scale factor 
at the time the strings formed. The extra term grows 
as cx 1/aH as the universe expands, ruling out a scaling 
solution. Eventually, when aH ~ e 2 T s a s , it starts to have 
a significant impact even on subhorizon scales. 

On the other hand, zero-temperature quantum vacuum 
fluctuations are still compatible with scaling. When their 
contribution is included, Eq. (|67p becomes 



G(k)/H 



\{H£)(k max /H) 3 /K, 



for k < k max /H, 
for k > k max /H. 

_ (69) 

where the coefficient of the extra term is automatically 
constant. This suggests that the vacuum fluctuations 
do not rule out scaling but modify the scaling solution. 
While this does not exclude the possibility of observable 
effects, it means that they would necessarily be quite 
subtle. 

The effects of the scaling violation (|55|) in the thermal 
case can be seen by considering the number of "infinite" 
strings. In the cosmic string terminology, this means 
strings that extend through the whole Hubble volume. 
The number of infinite strings at any time can be esti- 
mated by calculating the average string density in the 
Hubble volume, 



AitR 3 



\£\<R 



d 3 xp(x), 



(70) 



where R = l/H is the horizon size. Its typical value is 
given by the square root of the variance, 



avg/ 



/ \ 1/2 

-4-3 / d 3 xd 3 y2G(x-y)) 

3 / r°° f ]h 
— I / —(kRcoskR-smkR) 2 G(k) 



1/2 



dk k 2 G{k) 



1/2 



(71) 



where the contribution from k > 1/R = H has been 
ignored because it is the same as from the Kibble-Zurek 
mechanism and therefore small. Substituting Eq. ([58]) 
and ignoring numerical factors of order one, one finds 
that the average string number density in the observable 
universe is typically 



e 2 T s a s H 3 



(72) 
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This is due to strings that cross the horizon, because the 
net contribution from any loop that is fully inside the 
observable universe is zero. It be possible to construct a 
distribution of small loops that would reproduce Eq. | |72|) 
as a surface contribution, but Eq. (|7ip indicates that it 
arises from long- wavelength modes with k < 1/R. The 
most natural interpretation is therefore that it is given by 
infinite strings that extend through the whole observable 
universe. The typical number of infinite strings N Tms is 
then obtained by multiplying the density by the cross- 
sectional area irR 2 ~ 1/H 2 of the observable universe, 



iV rn 



Pirns /H 



e 2 T s a s 
aH 



(73) 



As the universe expands, this number grows. There 
would, therefore, be more and more infinite strings. 

Today, with a — 1 and H = Hq, the number of infinite 
strings would be 



N 



' e 2 T s a s 



(74) 



Had the relevant long-distance gauge field modes been in 
thermal equilibrium at the time of string formation, their 

1/3 

temperature would have been T s — T c w Tcmb/oss, 
where ~ 100 is the number of relativistic degrees of 
freedom at that time and Xcmb ~ 2.7K is the tempera- 
ture of the cosmic microwave background radiation. In 
this case the number of infinite strings today would be 
huge, iV ~ (IcMBZ-ffo) 1 ^ 2 ~ 10 14 , in massive conflict 
with observations. 

In inflationary cosmology, the long-wavelength super- 
horizon modes would not have been in thermal equilib- 
rium at the time of string formation, as was discussed 
in Section IV Bl Instead, their temperature would be 
T s = e~ N T[ n i, where N is the number of e-foldings and 
Tin; is the pre-inflationary temperature. Therefore we 
can write Eq. (|T4")) as 



N 



' e 2 T\ n \a 

1 1 



(75) 



where a; n i = e a s is the scale factor at the start of infla- 
tion. In order for inflation to solve the horizon problem, 
the comoving horizon size at the start of inflation must 
have been larger than today, i.e., Hi n ia ln i < Hq. If this 
bound is saturated, one finds 



N 



(76) 



In order to have inflation, the temperature Tj n i can not 
have been so high that is would have dominated the en- 
ergy density. This gives a bound 



H, 



> 9* 1 vox 



(77) 



where g* is the number of degrees of freedom and Mpi = 
(87rG) -1 / 2 is the reduced Planck mass. This bound is 
saturated if the temperature was initially even higher and 
inflation started when radiation became subdominant. In 
this case, one finds 



No 



3 2 M, 



PI 



9 



1/2 



T- 



(78) 



This can be as high as a dozen for GUT scale inflation, 
and in low-energy inflationary models it can be much 
higher. 

The presence of these infinite strings may affect on 
the evolutions of the string network and have indirect 
observable consequences. Their contribution to the total 
energy density would be 



JtNo_ 

3M 2 , : 



(79) 



which can be significant if Nq is large. Furthermore, they 
would single out a preferred direction, that of /? aV g in 
Eq. ([70]) . thereby breaking the isotropy of the universe. 
It may also be possible to observe them directly, and 
count them, using the Kaiser-Stebbins effect [Toj | . which 
creates discontinuities in the temperature of the cosmic 
background. 



VIII. CONCLUSIONS 

The results in this paper show that gauged cosmic 
string formation is dominated by trapping of gauge "mag- 
netic" flux at long distances, and the gauge field configu- 
ration becomes imprinted in the string network. Because 
of the preceding period of inflation and the masslessness 
of the gauge field, this can give rise to correlations over 
massively superhorizon scales at the time of string forma- 
tion and can therefore have potentially observable conse- 
quences. 

The detailed form of the string correlations was calcu- 
lated for thermal and vacuum initial states, correspond- 
ing to different cosmological scenarios. The correlations 
are stronger in the former case, but they are also present 
in the latter, cosmologically more realistic case, which 
describes a phase transition following a long period of 
inflation. In that case, these correlations are not gener- 
ated by inflation, and that the findings are therefore not 
in conflict with the well-known result [30| that inflation 
does not generate long-range magnetic fields. Instead, 
the correlations are a property of the vacuum quantum 
fluctuations, which are turned into physical string exci- 
tations by the phase transition. The role of inflation is 
simply to remove all other fluctuations, thermal or oth- 
erwise, which would lead to stronger correlations. The 
vacuum result is therefore not specific to inflationary cos- 
mology, and be seen as a general lower bound given by 
the uncertainty principle. 
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The results were confirmed by classical lattice field the- 
ory simulations. Although this is a widely used approach, 
its validity in the zero-temperature case is uncertain, be- 
cause the whole effect originates in vacuum fluctuations 
are is therefore fully quantum mechanical. Interestingly, 
the predictions can also be tested in superconductor ex- 
periments 44], which are obviously fully quantum me- 
chanical and therefore avoid this potential pitfall com- 
pletely. 

The calculations in this paper were carried out in the 
Abelian Higgs model, as is the case with essentially all 
studies of cosmic strings. The qualitatively new be- 
haviour is therefore not a consequence of adding new de- 
grees of freedom or modifying the microscopic theory in 
other ways. Instead, it arises because instead of impos- 
ing ad hoc initial conditions, the fields are set to be in 
their thermal equilibrium or quantum vacuum state. Be- 
cause the arguments used were based on general physical 
principles, the results should also generalise to essentially 
any model of local cosmic strings that carry a gauge flux, 
including cosmic superstrings formed in brane collisions. 

The observational effects of these correlations are not 
yet known. It was shown that sufficiently strong ther- 
mal fluctuations will violate the usual scaling assump- 
tion, and even at zero temperature, the precise scaling 
solution would have to change. This may well affect the 



cosmic string contribution to the CMB power spectrum. 
There may also be totally new effects. For instance, one 
can speculate whether the preferred direction of the infi- 
nite strings produced by flux trapping could be observed. 
Also, the difference between string networks in different 
Hubble volumes are not random, as in the conventional 
pictures, but correlated over infinite distances. Since 
they influence the expansions of the universe locally in 
each Hubble volume, they can give rise to potentially ob- 
servable density perturbations. Determining these signa- 
tures will require large-scale numerical simulations of the 
string network evolution. A key lesson from this paper is 
that the usual random phase initial conditions, which are 
a concrete realisation of the Kibble-Zurek mechanism, are 
insufficient and, instead, full quantum or thermal initial 
conditions should be used. 
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